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The 2-leg t—J ladder forms a spin liquid at half-filling which evolves to a Luther-Emery liquid 
upon doping. Our aim is to obtain a complete phase diagram for isotropic coupling (i.e. rungs and 
legs equal) as a function of electron density n and the ratio J/t. Two known limiting cases are: 
n < 1/2 which is a single band Luttinger liquid and small hole doping 8 <C 1 for J/t — > which 
is a Nagaoka ferromagnet. Using Lanczos techniques we examine the region between the Nagaoka 
and Luther-Emery phases for 1 > n > 1/2. We find evidences for gapless behavior in both spin and 
charge channels for J/t < 0.3 consistent with Luttinger liquids in both bonding and anti-bonding 
bands (i.e., C2S2). This proposal is based on the behavior of spin and charge correlation functions. 
For example the hole-hole correlation function which displays hole pairing at larger J/t, shows 
hole-hole repulsion in this region. As a further test, we examined the dependence of the energy on 
a relative phase shift between bonding and antibonding bands. For J/t < 0.3 this is very weak, 
indicating a lack of pairing between these channels. 

PACS numbers: 71.27.+a, 71.30. +h, 74.72.-h 



I. INTRODUCTION 

The surprising,discovery of high-T c materials by Bcd- 
norz and Mulleiu has sparked renewed interest in low- 
dimensional strongly correlated quantum systems. A 
striking feature of these materials is that they show sim- 
ple long range antiferromagnetic (AF) order at low tem- 
peratures when they are not doped with holes. 

Between the well-known one-dimensional systems and 
the difficult two dimensional system, coupled chains 
(or ladders) are interesting intermediate systems. The 
most striking feature of 2m-leg ladders is the appear- 
ance of a spin-gap at half-filling and small dopingptl 
For instance compounds such as Sr„_iCu„+i02n have 
been shown to be well described by a lattice of cou- 
pled f„-leg ladders, with v n = (n + 1)/2lE3. An- 
other example containing two-leg ladders is the system 
Sri4_ K Gaj r Cu2404iEJIlj which is a material with doped 
laddersriri and where sup^aapductivity (under high 
pressure) has been observedliThZl. 

An important step is the .determination of the phase 
diagram. Balents and Fisherlia have computed the phase 
diagram of the Hubbard-ladder in the weak coupling limit 
U — > in the frame of bosonization and renormalization 
group theory. To distinguish the different phases they in- 
troduced the notation CxSy for labeling phases with x 
gaplpss charge and y gapless spin excitations. Noack et 
al.t3 investigated then the different phases numerically 
with Density Matrix Renormalization Group (DMRG) 
methods for intermediate coupling and found good agree- 
ment with the analytical results of the weak coupling 
limit. The main feature is the existence of two regions, 
the Luther-Emery (LE) region with a spin gap (C1S0) 
and the Luttinger-Liquid (LL) with one gapless excita- 
tion (C1S1). pp 

The strong coupling t-J modelE3'Eil has also been stud- 



ied extensively by different authors. It is given by the 
Hamiltonian 

H = t 1 _Y J (4j,^2j,a + h.C.) (1) 
Ija 

3 ^ ' 

+ J II ( S/ < S ~ \ n U+i n ij ) > ( 2 ) 
lj v / 

where the index I £ {1,2} refers to legs and j to 
rungs. The c operators denote the fermion operators 
with projection onto the singly occupied states, i.e. 
Qj.o- = Qj,o-(l ~ nij.-o-)- In this paper, we report a de- 
tailed study of the phase diagram for isotropic coupling 
(i.e. rungs and legs equal) as a function of electron den- 
sity n and the ratio J/t. 

In the strong coupling limit Jj_ 3> J\\ , the spin-gap can 
be easily represented. At half-filling, the ground state is 
formed by spin-singlets lying on each rung. Turning over 
one spin gives a triplet on the corresponding rung. The 
energy difference between the two states, the spin-gap, 
is A ~ Jx . Numerically, it is found that the spin .gap 
in the isotropic case J|| = Ji reduces to A k Jx/2tTQ. 
The spin-gap of a doped ladder for J± ^> J|| is due to a 
qualitatively different process. Some singlets have been 
replaced by hole-pairs moving along the ladder with a 
renormalized hopping ti\ . The lowest excitation arises by 
breaking one hole pair (or equivalently a singlet bond) 
and putting each unbound spin on well-separated rungs. 
Turning over one of the separate spins gives the lowest 
triplet excitation with A i=s J± — 2t± — 2t\\ leading to 
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a finite spin-gap (LE). At the isotropic point J± = Jm 
and for J/t > 0.3 the spin-gap still exists where now the 
bound holes share neighboring rungs and stroua-sjnglet 
correlations are measured on the remaining sites&EE The 
spin-gap is thus a discontinuous function of doping in the 
vicinity of half-filling. For higher values of the parameter 
J/t, holes and electrons separate completely at all doping 
levels. 

Thus, the appearance of the spin gap in a doped ladder 
is directly correlated with the formation of hole pairs, and 
it is an interesting problem to study their stability. Poil- 
blanc et a/.cJ used a numerical method based on Lanc- 
zos algorithms, which has allowed them to distinguish 
between the gapped and gapless rtegions in the phase di- 
agram. Hayward and Poilblanccil determined the non- 
universal correlation exponents K p defining the behavior 
of the long-range correlations. They found at low electron 
density n < 1/2 the system is in a LL phase, while for a 
higher electron density a gapped phase with hole pairs is 
stabilized. A large region of this gapped phase exhibits 
dominant superconducting correlations. The boundary 
of these two phases was determined to be n ~ 1/4 where, 
in the band picture, the Fermi energy Ep just touches 
the antibonding band. 

However, some parts of the phase diagram are still un- 
clear. First, the finite-size scaling process to determine 
the spin-gap region does not give clear results for val- 
ues of J/t < 0.3. Second, the radius of the hole pairs 
increases with decreasing J/t, and it is not clear whether 
it diverges for a particular value, and if a transition to a 
gapless phase occurs for small J/t. For very small values 
of J/t, the t-J model is similar to the Hubbard model 
with U — ► oo. Since for two coupled chains, all spin con- 
figurations with a fixed S z are coupled by hole hopping 
as in a two-dimensional syistem, the essential condition 
for the Nagaoka theoremc3 is fulfilled. A ferromagnetic 
phase occurs for very small values of J/t and low hole 
doping. In this phase, no spin gap occurs and holes repel 
each other. With increasing J/t, the ground state rapidly 
evolves into a singlet state. 

This part of the phase diagram will be extensively stud- 
ied on the basis of exact diagonalization results for small 
clusters, typically 10-rung ladders. Finite-size effects can 
be important for such systems, and it will be tried to 
minimize them as much as possible. 



II. BOUNDARY CONDITIONS 

A Lanczos algorithm will be used to investigate the dif- 
ferent phases. With current computers it is possible to 
investigate 2 leg ladders of length L = 10 at any filling. 
In order to carry out a systematic analysis for different 
doping, the boundary conditions (BC) must be chosen 
carefully. 



A. Definition of Boundary Conditions 

Usually, the BC are defined in the non-interacting 
limit of the model where the Hamiltonian is exactly 
described by two parallel bands E±(k) also labeled 
E{k x ,k y ) with k y € {0,7r}. Since the system is finite, 
the k(= (k x , k y ))-v&hies belong to a discrete set. In gen- 
eral k x = ^-l + (j>, where a priori (j> needs not be the same 
for bonding (k y — 0) and antibonding (k y — n) branches. 
BC giving either fully occupied or empty single particle 
orbitals in the non-interacting case, called closed shell 
(CS) boundary conditions (CSBC), are chosen. 

If this condition is not fulfilled, the ground state is de- 
generate This favors the pairing of spins at Ep and en- 
hances pairing instabilities. This configuration is called 
open shell boundary conditions (OSBC). 

The definition of the BC for the bonding and antibond- 
ing operators 



bj„ = -/=( c lj> + C 2j>) ! (3) 

O-ja = -y=(cij,a-C2j,a), (4) 

are given by the set of equations (the spin index is droped 
for simplicity) 

TbtT- 1 = e<*&t +11 l<j<L, (5) 

TotT" 1 = e**o} +1> l<j<L, (6) 

Tb^T- 1 = e^bl, (7) 

Ta\T~ x = e i ^ +m ^a\. (8) 



For m = the usual BC with the geometry of a ring 
are recovered. They will be generally referred to as 
RBC(0i). (f>L = 0(71") are the most-used phases, which 
will be called periodic (antiperiodic) boundary condition 
or PBC(APBC). For m — 1 the operators at the end of 
the legs fulfill the relationship 

Tc^T- 1 = e^cU, (9) 

r4 i r- 1 = e i Mi, (10) 

corresponding to the geometry of a Moebius band. They 
will be called Moebius Boundary Conditions and denoted 
by MBC(^L). MBC(0) means that periodic boundary 
conditions for bonding states and anti-periodic boundary 
conditions for antibonding states are taken and vice- versa 
for MBC(7r). In real space they can be viewed as a way 
to prevent antiferromagnctic frustration in some cases as 
schematically shown in Fig. |l| for a 10-rung ladder with 
two holes. For a phase differences mir other than or w 
the translation of one particle leads to two states, each 
one with one particle on each leg. This case will not be 
considered further. 

By introducing these BC it is possible to get CSBC for 
all doping with an even number of holes. In the situation 
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where two CSBC are possible, the one giving the minimal 
energy in the t-J model will be generally preferredEi 
The Fourier transform 

h X'' kJ ./j- (ii) 

j 

with k = {k x , k v ) and j = (ji, J2) will be computed con- 
sistently by taking k x as a function of k y £ {0, 7r}, 

f *rj RBC, 
*»(*») = 1 , (12) 

where Z is an integer. The two possible sets of k- values 
are summarized in Fig. |^. 

III. NAGAOKA PHASE 

As discussed in the introduction, a Nagaoka phase!! 
for J/t — ► is expected at low doping. Numerically, 
TroyerB has shown that the ground state of L > 4-rungs 
t-J ladders with a doping of two holes is a saturated fer- 
romagnet at J/t — 0. In the following, the doping and 
J/t dependence of this phase will be investigated. 

The analysis of that phase is simplified by first con- 
sidering the subspace of the completely polarized states 
(5*°* = N s /2 with N s the number of spins) which is 
equivalent to a spinless fermion system. The eigenval- 
ues are exactly given by the band picture E±(k) = 
— 2tii cos(k x ) ± t± and the eigenstates are direct prod- 
ucts of bonding and antibonding states. For finite clus- 
ters CSBC can be sometimes obtained with different BC. 
For instance, for a 10-rung ladder with 2 holes and with 
tj_ = t\\ both APBC and MBC(O) allow a CS configura- 
tion with the same ground state energy as schematically 
plotted in the upper graphs of Fig. |3J. For higher J/t, 
when the ground state is a singlet, the corresponding 
CSBC configurations are shown in the lower graphs of 
Fig. H where black circle stand for fully occupied states. 
Both MBC(O) and MBC(tt) are possible. 

In Fig. [I| and [| the ground state energy is shown for dif- 
ferent cases. In Fig. [| the upper graph shows the lowest 
energies for a 5-rung ladder with one hole and the lower 
graph for a 10-rungs ladder with two holes. In Fig. || the 
lowest energies for a 5-rung ladder with two holes and for 
a 10-rung ladder with four holes is shown in the upper 
and lower graph, respectively. Each curve corresponds to 
different BC. For each graph, the transition to the satu- 
rated ferromagnet where the energy is J/t-independent is 
clearly seen for CSBC. The two other BC do not have CS 
and the corresponding energy states are singlets showing 
a linear J/t-dependent energy lying above the energy of 
the fully polarized state. This case will not be considered 
further. 

A crude finite size scaling can be made by extrapolat- 
ing the critical values of a 5-rung and a 10-rung ladderEJ. 
At a doping of 6 = 0.1, for L = 5 the lowest critical 



value is given by APBC with value (J/t)s ~ 0.052 while 
for L = 10 both APBC and MBC(0) gives the same crit- 
ical value of {J/t) 10 — 0.042. The infinite extrapolation 
gives {J/t)oo ~ 0.032. For 8 = 0.2, the two values are 
(J/t) 5 0.078 and (J/t) 10 ^ 0.056 respectively, and the 
infinite extrapolation gives (J/t)oo — 0.034. For higher 
doping, i.e. for a 5-rung ladder with three holes and for 
a 10-rung ladder with five and six holes with J/t > no 
ferromagnetic ground state has been found. 

In conclusion the diagonalization of small clusters 
shows signs of saturated ferromagnetism up to a doping 
of S ~ 0.2 and for small values of J/t {J/t ~ 0.033). 



A. For J/t > {J/t)c, Sign of Partial Ferromagnetism? 

For some clusters, the transition from the Nagaoka 
phase to the singlet phase does not occur immediately, 
but passes through different phases with an interme- 
diate spin < S < N e /2. The lowest energy of the 
rih = 2, 10-rung ladder for J/t < 0.09 is always given by 
APBC (with the exception of a very tiny region close to 
J/t = 0.04). In that region, the spin of the corresponding 
ground state passes through the finite values S — 5 for 
0.043 < J/t < 0.057 and S = 1 for 0.058 < J/t < 0.062. 
For J/t > 0.62 the ground state is a singlet with a fi- 
nite momentum of k = tt/5 and a singlet at zero mo- 
mentum for J/t > 0.083; they are mentioned at the 
bottom of Fig. || For a 10-rung ladder with = 4 
(Fig. ^|), some partial ferromagnetism also occur. For 
0.057 < J/t < 0.072 the spin is S = 2 with MBC(0) giv- 
ing the lowest energy. For J/t > 0.072, the lowest energy 
is given by PBC at k x — and is a singlet. 

The spin-spin correlations for the 10-rung ladder with 
n/i = 2 and for S z = S are plotted in Fig. |[ The corre- 
lations in real space are plotted in the uppermost graph 
where the sites on the first leg are labeled 1 < j < 10, and 
on the second leg 11 < j < 20 as pictured in Fig. |l|. For 
J/t > 0.05, they show an alternating behavior around a 
ferromagnetic value indicating antifcrromagnetic correla- 
tions. However, their Fourier transform, plotted in the 
lower graph, do not clearly indicate a continuous process. 
For J/t = 0.07 the correlation functions have a maximum 
in the branch k y = tt indicating that the sum of the inter- 
band scattering processes is greater than the intraband 
processes and thus, that the correlations between both 
bands are important. The correlations at other values 
of J/t display a maximum in the k y — branch showing 
that the intraband processes are now favored. 

These effects are not found for small systems (L=5). 
It is thus tempting to conclude the existence of narrow 
region with partial ferromagnetism phase in the thermo- 
dynamic limit of the n-J phase diagram. However the 
present results are not enough to draw a clear conclusion 
on this subject and the existence of these phases is still 
open. 
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IV. EXISTENCE OF A C2S2 REGION? 

The spin gap of a doped system is intrinsically bound 
to the formation of hole pairs. Moreover, it is known that 
for very low J/t, when the system is ferromagnetic, holes 
repel each other. The question of whether the transition 
between repulsive and attractive holes occurs at the fer- 
romagnetic to paramagnetic transition will be addressed 
in this section. Numerical evidence shows that a region of 
holes with repulsive residual interactions, and thus with 
gapless spin excitations, occurs between the ferromag- 
netic and the LE phase. 



A. Hole— Hole Correlations 

A first insight into this question can be ob- 
tained by looking at the hole-hole correlation function 
(nh(l)rih(j)). In Fig. ^ the correlations are plotted 
for different values of J/t and different fillings, namely 
nh = 2 with MBC(tt), n h = 4-fvith PBC, n h = 6 with 
MBC(O) and n h = 8 with PBCEL For n h = 2 and n h = 4 
the correlation for the ferromagnetic state computed with 
APBC is also shown (circles). The uppermost graphs in 
the figure show the correlation functions in real space 
with the site convention pictured in Fig. [j]. Each corre- 
lation is normalized such that (rih(l)nh(l)) = 1. 

The lower graphs show the Fourier transform 

(n h (k)n h (~k)) oc Af(k) = £ e* kj (n h (l)n h (j)), (13) 

j 

where k is defined in Sec. |l[ The normalization cor- 
responds to jV(0) = rih (out of the graphs). In Fourier 
space, the points in the k y = 0(tt) branch are given by 
the left(right) curves in each graph. 

In real space, each graph shows the same behavior. For 
instance for n/j = 2(6 = 0.1) at J/t — 0.1 the second hole 
is found to sit at the farthest point from the first hole. 
When J/t is increased the hole-hole correlation for hav- 
ing two holes on the same rung increases. When J/t = 1 
the hole-hole correlation clearly shows their bound char- 
acter. 

In the large Jj_,t_i_ limit, with J± > 2t j_, the hop- 
ping of hole pairs between rungs is given in second or- 
der perturbation theory by t|| = —2t?,/(J± — 2t±). As 
the hole pairs act as hard-core bosons, they repel each 
other in order to gain the maximal kinetic energy. This 
can be viewed in the lower graphs plotting the Fourier- 
transform, where peak in the k y — branch appear. This 
feature clearly appears in the four different graphs. At 



the same time, the total weight in the k y = it decreases. 

This behavior raises the question of the existence of an 
intermediate phase below the LE phase. In fact, for two 
holes and values of 0.05 < J/t < 0.2 the system seems to 
be in a phase where holes repel each other. According 
to the above picture, it would imply spin-gapless excita- 
tions. However, it could also be a pure finite size effect in 



that the radius of a two-hole bound state is greater than 
the length of the ladder sample. Thus, finite size effects 
strongly complicate the interpretation. The above cor- 
relations are not enough to conclude definitively to the 
existence of a new phase. 



B. Spin— Spin Correlations 

The spin-spin correlation of the systems are plotted in 
Fig. M. The uppermost graphs show the spin-spin corre- 
lations (S z (l)S z (j)) in real space. They clearly show the 
antiferromagnetic ordering of the spins along the chain. 
For rih = 4 and — 6 and for low value of J/t, the cor- 
relation across the rung (S Z (1)S Z (11)) is ferromagnetic 
not antiferromagnetic. Increasing J/t however stabilizes 
the system to have antiferromagnetic ordering across the 
rung. 

This feature is also clearly emphasized by the Fourier 
transform, plotted in the lower graph of Fig |8|, which 
have their maxima in the k y = branch. Their Fourier 
transforms A4(k) also gives information on the different 
scattering processes occurring between the different spins 
in the band picture. In fact, they are proportional to the 
equal-time correlation function, 



M(k) oc (S z (k)S z (-k)) = ^|(a|5 z (A;)|0)| 2 



(14) 



where the sum is taken over all eigenstates of H 
and where |0) denotes the ground state. For the 
case 7ih = 2, the MBC(7r) has bonding states filled up 
to /cp.b = \ and antibonding to fcp.a = \ as schemat- 
ically shown in Fig. ^. The lowest energy in- 
traband excitations in the bonding and antibonding 
band occur at k 6 = (2/c F ,b + x> °) = (tF> °) and at 
k a = (2fc F , a + 22L,0) = (¥,0), respectively. The corre- 
sponding curves in Fig. show a maximum at kf, while 
no peak appear at k a showing that these processes do 
not dominate in the sum of Eq. [l4|. Moreover, these 
processes are quite small compared to the interband pro- 
cesses plotted in the right part of the graph. According 
to the band picture the lowest energy interband processes 
are at k a & = (&F,b + ^F,a + \ = ff , 7r). The curve display 
a maximum at that point dominating all the other pro- 
cesses and especially the intraband scattering ones. It 
leads to the interpretation that the bonding and anti- 
bonding particles are strongly correlated at Ef- In the 
next section a picture will be proposed in which the gap- 
less phase is a manifestation of the absence of correlation 
between bands. In such a picture, the correlation for 
nh = 2 is characteristic for a gapped phase. 

For nh = 4(PBC) the band picture in the non interact- 
ing limit predicts dominant intra-bonding band scatter- 
ing processes at kj, = (tt, 0) and intra-antibonding band 
scattering processes at k a = (3tt/5, 0). Dominant in- 
terband scattering processes occur at k a f, = (4tt/5,tt). 
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These values correspond to the (local)maxima of the cor- 
relation in agreement with the simple non-interacting 
band picture. 

For rth = 6(MBC(0)) the dominant scattering pro- 
cesses are expected at k& = (7r, 0), k a = (2tt/5, 0), 
and k Qb = (77r/10,7r). For n h = 8(PBC) the domi- 
nant scattering processes are expected at kf, = (-7r, 0), 
k a = (7r/5,0), and at k a b = (37t/5,7t). The correspond- 
ing graphs display corresponding maxima in good agree- 
ment with the band-picture predictions. 

In Fig. ^ the maxima of the correlation functions for 
the k y = (circle) and the k y = it (square) branches are 
plotted. For — 4, 6 and Uh = 8 they show a crossover 
from a region with dominant intraband scattering to a re- 
gion with dominant interband scattering. This suggests 
that the system has uncorrelated bonding and antibond- 
ing bands at Ep such that the low energy physics is anal- 
ogous to that of two one-dimensional systems, with spin 
and charge gapless excitations (C2S2). A rough criteria 
for the transition can be defined by taking the critical 
J/t at the crossover. This yields (J/t) c ~ 0.32,0.34,0.26 
for 5 — 0.2,0.3, and 5 — 0.4, respectively. For nj, = 2 
the interband scattering is always much bigger than the 
intraband scattering yielding to the conclusion that no 
C2S2 phase occurs for that particular filling. This em- 
phasizes that care must be taken in interpreting the sim- 
ple hole hole correlations. 

C. (Anti-) Bonding Pair Correlations 

Defining the singlet 5? and triplet f} , f} a , a e {1,1} 
creation operator on the rung i with 

= "^( C iiT C 2U C HJ. C 2it)' (^) 
Ti a = c \ia c 2iai (1^) 

the different pair correlations between bonding and anti- 
bonding states inr-the subspace of single occupied states 
can be written asEj 

(blA&Wl) = (18) 

In the C2S2 phase, when the bonding and antibond- 
ing states are not correlated, all pair correlations are 
shortrange. In the spin-gapped region, the ground state 
is characterized by the formation of interchain singlet. 
Thus, long range singlet-singlet (SS) pair correlations 
will appear while the triplet-triplet (TT) pair correla- 
tions will remain short-range. 

For an undoped system, the pair correlations are finite 
only on the same rung j = i, where they are proportional 



to the number operator for singlet or triplet states, and 
depend only on the ratio J±/J\\. For J|| = 0, the singlet- 
singlet pair correlation at j = i is 1 while the triplet- 
triplet pair correlation at j = i is 0. At the isotropic 
point, the values for the undoped 10-rung ladder are 

(£t(i)£(i))~ .7, (19) 
(ft(l)Tb(l)) = (ft(l)f T (l))~0.1. (20) 

In a doped system, the pair correlations at the isotropic 
point depend on the ratio J/t. In Fig. [k] (5^(1)5(1)) 
and (f t (l)f (l)) are plotted as function of J/t. The 
case n h = 2 shows that (5J (1)5(1)) > (f£(l)f (l)) for 
all plotted values J/t. It favors the scenario that a 
10-rung ladder with two holes has a gap for very low 
J/i-values in agreement with the discussion of the spin- 
spin correlation. With increasing J/t the singlet num- 
ber (5t (1)5(1)) rapidly increases while the triplet num- 
ber (Tq (l)Tb(l)) decreases. On contrary, for = 4(6) a 
cross-over between the singlet and triplet number occurs 
at J/t ~ 0.23(0.17). This is in agreement with the spec- 
ulated gapless region at low J/t. The case with = 8, 
exhibits a favored rung singlet configuration. 

In Fig. [ll] the pair correlation at different rungs 
is shown as a function of J/t for the different 
cases. To make a consistent comparison of the 
long-range order for various J/t the pair correla- 
tions are normalized according to the singlet number, 
(S](l)S(j)) r = (5t(l)5(j))/(5t(l)5(l)). For n h = 2 the 
pair correlations uniformly tend to zero when J/t de- 
creases to zero. The solid and dotted lines show the 
derivative (forward difference) of the curves for j — 2 and 
j = 3. Both show a peak around J/t = 0.3 indicating 
the value of J/t below which the pair correlation become 
short-ranged. The curve does not show any different be- 
havior for J/t > 0.1. Below this value, it has been seen 
that the system has a transition to a ferromagnetic phase 
at J/t = 0.032. Between these two points some sign of 
partial ferromagnetism has been discussed but no direct 
evidences for a C2S2 phases has been observed. 

The striking feature of the graphs for — 4, 6 and 
rih = 8 is the qualitatively different behavior of the 
different curves at small J/t. The normalized short 
range pair correlation (5 t (l)5(2)) r shows a minimum 
at J/t = 0.30, 0.45, and 0.50 for n h = 4, 6, and 8, re- 
spectively. The solid and dashed/dotted lines show the 
derivative of the different curves. Maxima in the deriva- 
tives appear in both cases for 0.25 < J/t < 0.5 below 
which value the correlation get short range. 

A critical value will be defined to be at the point 
where the long-range pair correlation (5 t (l)5(5)) J , 
change the slope of its curve. This gives the values 
(jft) c ~ 0,20,0.25,0.25 for 5 = 0.2,0.3, and 0.4, respec- 
tively. 
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V. TRANSFORMED HAMILTONIAN 

In this section, the physical picture for the C2S2 
phase based on the absence of phase correlation between 
the bonding and anti-bonding bands is examined. For 
J/t ~ 0.5 phase coherence between the bands at Ep ex- 
ists leading to a LE liquid with the spingap. The ques- 
tion whether this phase coherence disappears for small 
J/t leading to two independent LL (C2S2) will be inves- 
tigated^. The terms in the Hamiltonian coupling both 
bands are of the form b\ x .% i<7 /afc 3 ,<r'Ofc 4 ,o- + b.c. They 
must be irrelevant near Ep in the gapless LL phase and 
become relevant when J/t increases. To investigate this 
we transform the Hamiltonian by introducing a relative 
phase between the bonding and antibonding operators, 
so that only terms of the above form are affected, then 
the dependence of the ground state on this phase differ- 
ence will be studied. 



A. Relative Phase Transformation 

The transformation is defined by 

H = PH(b, a)P -► PH{e l4 ">b, e i0 "a)P = H t , (21) 

where P is the projector onto singly occupied states and 
-ff(a, b) is the t—J Hamiltonian (without projection) writ- 
ten in the bonding and antibonding basis. The hopping 
term of the Hamiltonian is diagonal in the bonding and 
antibonding basis and thus does not change. The intro- 
duction of a phase shift will only affect terms of the form 

which occur in the Hj^ and H j ± part. They couple the 
bonding and the antibonding operators and are thus re- 
sponsible for the appearance of the spin gap. 

With 7 = (f>b — (j) a , The total transformation is sum- 
marized through the definition, 

PH( Clj )P — ► PH(c m )P = flt( 7 ), (22) 

where the transformed Hamiltonian H t (j) has a peri- 
odicity of 7r in 7. The site operators are transformed 
according to the rule c/j — ► cy,t with 

= \ ((e* + e^) Clj + (e**» - e^)c 2j ) (23) 

c 2j ,t = \ ((e* - e^)c l3 + (e**- + e^)c 2j ) . (24) 

The final transformed Hamiltonian can be written with 
three distinct terms 

where the first part i?t,h is the hopping part being the 
same as in the original Hamiltonian. The second term 



•fft,Jii involving magnetic coupling along the chain con- 
tains four distinct terms 

-fft,J|| = -^t,J||,l + #t,J||,2 + ft, J|| ,3 + #t,J||,4 (25) 

which are listed and discussed in Appendix [a|. The in- 
terleg magnetic coupling H tt j ± vanishes at 7 = 7r (see 
Appendix |A|) in agreement with the suppression of corre- 
lation between the bonding and antibonding operators. 



B. Ground State Energy 

The ground state may not be completely decoupled 
into a bonding and an antibonding part. Only the parti- 
cles with a k- value near Ep must be uncorrelated for two 
independent LL. By calculating the ground state energy 
of H t the influence of the phase shift on all bonding and 
antibonding operators is included and it is thus possible 
that the influence of the phaseshift is qualitatively sim- 
ilar for different values of J/t. In fact, in Fig. |l2] the 
ground state energy is plotted for a 10-rung ladder with 
4 holes and using PBC. By switching on 7, the energy 
increases for all plotted values of J/t. However, no clear 
qualitative change in the evolution of the energy vs 7 
for different J/t can be observed. This is a consequence 
of the fact that all particles (also away from Ep) are 
involved. It is interesting to consider the different corre- 
lation functions to look for signs if the system undergoes 
a phase-transition when 7 is altered. 



C. Correlation Functions 

In Fig. [ll| the hole-hole correlations for a 10-rung lad- 
der with 4 holes are plotted. The correlations are shown 
at different J/t for different values of 7. The upper 
graphs show the correlations in real space, and the lower 
their Fourier transform. For J/t < 0.3 only small changes 
in the shape of the correlations are observed when 7 in- 
creases. For J/t = 0.4 and J/t = 1.0 the shape of the 
correlation is clearly changed indicating state with re- 
pulsive interactions at 7 = 0.57T. The influence of 7 
for J/t = 1.0 is the most important. The rung corre- 
lation (nh(l)^h(H)) is clearly suppressed in agreement 
with the breaking of hole-pairs. The diagonal correlation 
("h(l)"h(12)) shows a peak at 7 = 0.27T which is sup- 
pressed for 7 = 0.5-7T. The peak in the branch k y = of 
the Fourier transform, characteristic of the homogeneity 
of the hole-pair repartition is also suppressed. 

In Fig. m the spin-spin correlation for the same clus- 
ter is plotted. The measured antiferromagnetism is quite 
sensitive to the phase shift. In real space the antiferro- 
magnetic correlation across the rungs is turned into leg 
independent behavior (rung ferromagnetism) . However, 
.M(k) measures all spins in the band pictures, therefore 
not only those at Ep, and a change in the correlation is 



G 



also compatible with the gapless phase. The most strik- 
ing picture of the spin-spin correlations is given by their 
Fourier transforms. As discussed above, dominant peaks 
in the bonding channel (k y — 0) are expected for the gap- 
less phase, while dominant interband (k y — tt) scattering 
processes are expected in the gapped phase. This behav- 
ior appears clearly in all graphs. 

In Fig. 15 the maxima of the spin-spin correlations as 
a function of 7 are shown for the different J/t values. 
They are at k = (it, 0) (circle) and k= ( ^ , tt) (square) 
(with the exception of J/t — 0.1, where the maximum 
is at k = (^p,0) when 7 > 0.4-7r), emphasizing the intra- 
band processes in the antibonding band. The crossing of 
the maxima for J/t > 0.4 is related to the destruction of 
the hole pairs. It occurs at 7 ~ 0.1 l7r for J/t = 0.4 while 
for J/t = 1.0 the crossover occurs at a higher 7 ~ 0.21-7T 
since the hole pairs are more strongly bound. 

The influence of the relative phase 7 on the singlet 
and triplet pair correlation functions is plotted in Fig. p"6] 
and |l7|. In the upper graph of Fig. [l6]the pair correlation 
(S't(l)S(j)) and (f£(l)f (j)) at J/t = 0.1 are plotted for 
different 7's. Pair correlations along the rung decrease 
rapidly to zero. At j = 1 the pair correlation measures 
the number of singlets (triplets) on the rung. The num- 
ber of singlets is much lower than the triplet number of 
the first rung, in agreement with the destruction of the 
singlet liquid state. By switching on 7, the singlet num- 
ber decreases while the triplet number increases. The 
inset of the graph shows the number operators as a func- 
tion of 7. In the lower graph, the singlet and triplet 
pair correlations are measured in the spin-gapped phase 
at J/t = 1.0. The singlet pair correlation starts with 
a higher value than the triplet pair correlation. It de- 
creases rapidly to a lower value and shows a finite value 
at j = 6. When 7 is switched on, the on-site singlet is 
suppressed while the triplet pair correlation is increased. 
The crossover occurs at 7 ~ 0.257T. This on-site behavior 
emphasizes the picture of destroyed singlets on the rungs. 

The long-range behavior is shown in Fig. [It]. where 
the normalized SS(TT) pair correlation at site j ' = 6 are 
plotted as functions of 7. For J/t — 0.1 the pair correla- 
tions remain nearly unaffected by the introduction of 7 
while for J/t = 1.0 the pair correlations are suppressed 
and reach nearly the same value at 7 — O.bir. 



VI. CONCLUSIONS 



such region has been seen. However, a very tiny region 
between the ferromagnetic and the LE phase may occur. 

For larger doping clear signs of a C2S2 phase have 
been observed. Two different sets of critical values were 
found. Those obtained from the study of the crossover of 
the maximal value of the spin-spin correlation function 
M(k) in Sec. 



IV B are consistent with the values from the 



study of the long-range singlet-singlet pair correlation of 
Sec. |IVC[ Their mean values are ( J/t) c ~ 0.26, 0.29, and 
0.25 for 8 = 0.2, 0.3, and 0.4, respectively. In Fig. ||, 
they are included in the previous., phase diagram pro- 
posed by Hayward and Poilblanco. In the very small 
J/t region, the Nagaoka phase is shown. At quarter 
filling, in the non-interacting limit Ep just touches the 
anti-bonding band while Umklapp processes occur in the 
bonding band. These features are difficult to simulate in 
finite clusters and will not be discussed here. 
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APPENDIX A: TRANSFORMED HAMILTONIAN 

The introduction of the relative phase between bond- 
ing and antibonding operators does not transform the 
hopping term Hh but the magnetic part is considerably 
modified. The first term of the magnetic coupling along 
the chain Equ. |25|, 



t.Jii.i 



= ^„E 



1- j(l-cos(2 7 )) J (S.uSi,.,, 



+S 2 iS 2 i+i - -(n li n li+1 + n 2i n 2l+1 ) 



is similar to the original Hamiltonian with a renormalized 
coupling constant which decreases as 7 € [0, n] increases. 
For 7 = 7r the coupling strength is half that of 7 = 0. 
The second term 



In this paper, an extensive study of the small J/t- 
region at low-doping using exact diagonalization results 
of small cluster has been performed. First the hole- 
hole correlation showed hole repulsion for small J/t for 
a certain doping range indicating the presence of a C2S2 
phase23. However, finite size effect strongly affect the re- 
sult and different correlations were introduced in order 
to examine other properties of the system. For a very 
small doping (8 < 0.1), no clear sign of the existence of 



Hi 



W11 ,2 



i(l-cos(27))(S M S w+1 - 



S 2 iSii+i - -(n u n 2l+ i + n 2i n ll+1 



is similar to the first part but with antiferromagnetic cou- 
pling along diagonals. 
The third part, 



7 



J N ,3 



= J\\ 



-(l-cos(2 7 )) 



-7( C L lc 2i ~ 4i lc l») (cL+1 1c 2»+1 - 4i+l lc H+l) 



features spin-dependent hopping terms involving a simul- 
taneous hopping on two neighboring rungs. 
The last part, 



H 



t,J|,,4 



-(sin(2 7 )) x 



[(s H - s 2 i)(4 i+1 ^c 2i+ i - 4i+i^ c H+i) 

~(n u - n 2 i)(c\ l+1 lc 2l+ i - c\ i+1 lc u+1 ) 

-\( C li lc 2i - 4 4 lc H)( n li+l - n 2i+l) 



is a combination of spin and hopping operators involving 
a spin dependent hopping on one rung if the neighboring 
rung is occupied. 

The interchain magnetic coupling after some algebra 
can be written as 



H 



t,jj 



J± E ( 1 - V 1 - cos ( 2 7))) f SiiSai - \m 



vanishing for 7 = n. 
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FIG. 1. A 10- rung ladder with two holes on the same rung. Due to MBC The hole pair when sitting on the same rung do 
not frustrate the Neel ordering. The labeling convention used in future graphs is shown. 
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FIG. 2. The two different sets of k values for RBC and MBC for a 10-rung ladder. For MBC, the set of k values at k y = tv 
is shifted by tt/L respectively to the k v = branch. 





FIG. 3. The CSBC for a 10-rung ladder with two holes. The upper graphs represents the filling of the band with APBC and 
MBC(O) for a saturated ferromagnet. The lower graphs shows the filling in the non-interacting limit for S = for MBC(7r) 
and MBC(O). Filled circles represent the doubly occupied states while white circles stand for empty states. 
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FIG. 4. The upper(lower) graph shows the data for a 5(10)-rung ladder with 1(2) hole(s) for different BC. The insets show 
the magnified region where the saturated ferromagnet appears. The total spin of the ground state (using APBC) for the 10-rung 
ladder is indicated at the bottom of the lower graphs. 
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FIG. 5. The upper(lower) graphs show the data for a 5(10)-rung ladder with 2(4) holes for different BC. The insets show 
the magnified region where the saturated ferromagnet appears. The total spin of the ground state for the 10-rung ladder is 
indicated at the bottom of the lower graphs. 
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FIG. 7. The instantaneous hole-hole correlations for a 10-rung isotropic ladder at different doping. The upper graphs give 
the correlation in the site representation while the lower graphs show their respective Fourier transforms. The convention of 
the labeling is that of Fig. [j] and ^. 
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FIG. 8. The instantaneous spin-spin correlations for a 10-rung isotropic ladder at different doping. The upper graphs give 
the correlations in the site representation while the lower graphs show their respective Fourier transforms. The convention of 
the labeling is the one of Fig. hi and H. 
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FIG. 10. The pair correlations (5' t (l)S'(l)}, (Tq (l)Tb(l)) as function of J/t for a 10-rung isotropic ladder at different doping. 
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FIG. 13. The instantaneous hole-hole correlation of a 10-rung isotropic ladder with 4 holes for different values of J/t and 
7. The upper graphs show the correlation in the site representation while the lower graphs show their Fourier transform. The 
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FIG. 16. The pair correlations (,V (l)S'(j)} and (T (l)To(j)) at different 7 values as function of the rung index for a 10-rung 
ladder with n h = 4 holes. The upper graph shows data for J/t = 0.1 and the lower graph for J/t = 1.0. The insets show the 
on-site pair correlations (S' t (l)S'(l)} and (Tjj (l)Tb(l)) as function of 7. 
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FIG. 18. The speculative phase diagram taken from the work of Poilblanc et al. with the new phases added at low value of 
J/t and n > 1/2. When n < 1/2 the Fermi energy crosses only the bonding band and a single LL occurs. 
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